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ON THE DYNAMICS OF A DEGENERATE PARABOLIC EQUATION: 
GLOBAL BIFURCATION OF STATIONARY STATES AND CONVERGENCE 

NIKOS I. KARACHALIOS t AND NIKOS B. ZOGRAPHOPOULOS* 


Abstract. We study the dynamics of a degenerate parabolic equation with a variable, generally 
non-smooth diffusion coefficient, which may vanish at some points or be unbounded. We show 
the existence of a global branch of nonnegative stationary states, covering both the cases of 
a bounded and an unbounded domain. The global bifurcation of stationary states, implies-in 
conjuction with the definition of a gradient dynamical system in the natural phase space-that at 
least in the case of a bounded domain, any solution with nonnegative initial data tends to the 
trivial or the nonnegative equilibrium. Applications of the global bifurcation result to general 
degenerate semilinear as well as to quasilinear elliptic equations, are also discussed. 


1. Introduction 

The mathematical modelling of various physical processes, where spatial heterogeneity has a 
primary role, has usually as a result, the derivation of nonlinear evolution equations with variable 
diffusion, or dispersion. Applications are ranging from physics to biology. To name but a few, 
equations of such a type have been successfully applied to the heat propagation in heterogeneous 
materials E3EUEUEU, the study of transport of electron temperature in a confined plasma 
129 . the propagation of varying amplitude waves in a nonlinear medium [717) (and El for linear 
Schrodinger equation), to the study of electromagnetic phenomena in nonhomogeneous supercon¬ 
ductors cmnninniETO and the dynamics of Josephson junctions unim, to epidemiology and the 
growth and control of brain tumors M- 

In this work we continue the study, initiated in E3 of the qualitative behavior of solutions 
of some degenerate evolution equations (involving degenerate coefficients). Work I55[ concerns 
the asymptotic behavior of solutions, of a complex evolution equation of Ginzburg-Landau type. 
Here we study the following semilinear parabolic equation with variable, nonnegative diffusion 
coefficient, defined on an arbitrary domain (bounded or unbounded) H C ]& N , N > 2, 

d t 4> — div(er(;r) V</>) — A </> + |^>| 27 </> = 0, ieSl, t > 0, 

(l.i) o) = 4>o(x), sen, 

4>\dn = 0, t > 0. 

Equation o can be derived as a simple model for neutron diffusion (feedback control of nuclear 
reactor) E3E9- In this case </> (which must be nonnegative) and er stand for the neutron flux and 
neutron diffusion respectively. 

The degeneracy of problem CEU is considered in the sense that the measurable, nonnegative 
diffusion coefficient cr, is allowed to have at most a finite number of (essential) zeroes, at some points 
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or even to be unbounded. The point of departure for the consideration of suitable assumptions 
on the diffusion coefficient is the work 1X21 , where the degenerate elliptic problem is studied: we 
assume that the function a : 9 —+ K satisfies the following assumptions 

(H a ) cr e L 1 1 oc (fl) and for some a £ (0,2), lim inf-,,^ \x — z\~ a a(x) > 0, for every z € 9, when 
the domain f 1 is bounded, 

(H™) (J satisfies condition (H a ) and liminf| a .|_ >0o \x\~ ,3 <r(x) > 0, for some (3 > 2, when the domain 
9 is unbounded. 

The assumptions ( H a ) and (' H c p) imply (see P2J Lemma 2.2]) that (i) the set of zeroes is finite, 
(ii) the function cr could be non smooth (cannot be of class C 2 , if a £ (0,2) and it cannot have 
bounded derivatives if a € (0,1)). Moreover, in the unbounded domain case the function cr has 
to be unbounded. The approach in CHI, was based on Caffarelli-Kohn-Nirenberg type inequalities 
(see O). For some recent results concerning these inequalities and their applications to the 
study of elliptic equations, we refer to PEnUHH! 

The physical motivation of the assumption ( H a ), is related to the modelling of reaction diffusion 
processes in composite materials, occupying a bounded domain 9, which at some points they behave 
as perfect insulators. Following G3 pg. 79], when at some points the medium is perfectly insulating, 
it is natural to assume that cr(a;) vanishes at these points. On the other hand, when condition 
(H’p 1 ) is satisfied, it follows from PS Lemma 2.2], that in addition, the diffusion coefficient has to 
be unbounded. Physically, this situation corresponds to a nonhomogeneous medium, occupying the 
unbounded domain f 2, which behaves as a perfect conductor in 9\B/j(0) (see [23 pg- 79] ), and as a 
perfect insulator in a finite number of points in B/{(0). Note that when <99 ^ 0, the function a(x ), 
need not be locally bounded. These conditions arise in various simple transport models of electron 
temperature in a confined plasma. See m for a discussion concerning the one-dimensional case: 
the electron thermal diffusion is density dependent such that it vanishes with density, rendering 
the problem singular. Note that in various diffusion processes, the equations involve diffusion 
c t(x ) ~ |a;|“, a < N: We refer to [ 31 ! EH f° r equations describing heat propagation. 

The main purpose of this work is to combine basic results from the theory of infinite dimensional 
dynamical systems and bifurcation theory, to give a description of the dynamics of (II.111 . We 
remark here the crucial role of the conditions ( H a ) and ( H “) on the “degeneracy exponents” 
a, (3 which give rise to necessary compactness properties of various linear and nonlinear operators 
associated to the study of O and its related stationary problem (a degenerate elliptic equation). 
We are restricted in the case N > 2 since the case N = 1, despite its similarities with the higher 
dimensional case with respect to the definition and properties of the appropriate functional setting, 
recovers also important differences. For the definition and properties of the related function spaces 
and detailed discussions on one dimensional versions of generalized Hardy and Caffarelli-Kohn- 
Nirenberg inequalities, we refer to mm- 

More precisely, the first part of the present work is devoted to some results concerning the 
existence of a global attractor. While the result in E3, for the complex evolution equation, 
concerns the existence of a global attractor in L 2 (9), here it is verified that the dynamical system 
associated to JEU is a gradient system, and that there exists a connected global attractor in the 
weighted Sobolev space Dq’ 2 (9; a), the closure of Cq°(Q) with respect to the norm ||0|| 2 li2 = 

f n a(x) |V</>| 2 . This space appears to be the natural energy space for 11.111 . The main result of 
Section 3, can be stated by the following theorem. 
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Theorem 1.1. Let 57 C R w , N > 2, be an arbitrary domain (bounded or unbounded). Assume 
that a satisfies condition (H a ) or (TTff) and 


0 < 7 < 


2 — a 


2{N - 2 + a ) 


:= 7 


Equation m defines a semiflow 


S(t) : Dl’ 2 (n,a) ^ Dl’ 2 (n,a), 

which possesses a global attractor A in -Dq’ 2 (12; a). Let £ denote the (bounded) set of equilibrium 
points of S{t). For each positive orbit <f> lying in A the limit set u>{(f>) is a connected subset of £ on 
which J : Z7q’ 2 (57; a) —> R, 

■= \ er(x)\S7(/>\ 2 dx \<t>\ 2dx + 27 1 + 9 \<t>\ 21+2 dx, 

the Lyapunov functional associated to S(t), is constant. If £ is totally disconnected (in particular 
if £ is countable), the limit 

z+ = lim <f>(t), 

t—►+oo 

exists and is an equilibrium point. Furthermore, any solution of m. tends to an equilibrium 
point as t —> +oo. 


Further analysis is carried out, regarding the bifurcation of the corresponding steady states 
with respect to the parameter A £ R. More precisely, we prove the existence of a global branch of 
nonnegative solutions for the equation 

—divloladVw) = A u — litl 27 it, in Q, 

( 1 . 2 ) 

u|an = 0, 

bifurcating from the trivial solution at (Ai, 0), where Ai, is the positive principal eigenvalue of the 
corresponding linear problem 

—divlcrladVrt) = Xu, in fl, 

(1.3) 

w|an = 0. 

This is the main result of Section 4, described by the following Theorem. 

Theorem 1.2. Let 57 C WL N , N > 2, be an arbitrary domain (bounded or unbounded). Assume 
that a satisfies condition (H a ) or (Hff), and 

n 2 -a 

<7< 2{N -2 +a)' 

Then, the principal eigenvalue Ai of UJS is a bifurcating point of the problem ED and C\ x is 
a global branch of nonnegative solutions, which ”bends’ ’ to the right of Ai. For any fixed A > Ai 
these solutions are unique. 


The technique leading to the global bifurcation result, is included in the general strategy of 
the approximation of solutions of a degenerate partial differential equation, by constructing an 
approximate sequence of solutions of nondegenerate problems. The approximation procedure has 
been successfully applied to evolution ED El, and to stationary problems E3I13ED and in the 
context of bifurcation theory EHHE2EH] 

One of the main difficulties arising, on the attempt to establish the global character of the 
branch of nonnegative solutions for ED, is that Harnack-type Inequalities are not valid in general 
(see [221 Remark 3.2]). This is a common fact for non-uniformly elliptic equations '.']2 . However, 
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we refer to |T] and the references therein for generalized Harnack-type inequalities, applied to 
degenerate elliptic equations. Distinguishing between the bounded and the unbounded domain 
case, we consider different families of approximate boundary value problems. 

When 12 is bounded, na Lemma 2.2 and Remark 2.3] implies that under assumption (Tt a ), 
there exists a finite set Z := {z \,..., Zk} C 12 and r,S >0, such that the balls of center z l and 
radius r, B r (zi), i = 1 ,,k, are pairwise disjoint and 

., (i) a(x) > 5\x — Zi\ a for x G B r (zi) n 12, i = 1,..., k, 

(ii) tr( x) > <5, for x G 12 \ (J i B r (zi). 

Moreover if a satisfies (' H a ), then er > 0 in 0, the set of zeroes of er Z a := {z £ ft : cr(z) = 0} is 
finite, and Z a C z ([IE]. Remark 2.3]). It is not a loss of a generality to assume that Z n = Z. 

For convenience and simplicity, in the bounded domain case, we consider as a model for the 
diffusion coefficient, the function 


(1.5) 


a(x) = |a;| a , a G (0,2), 


satisfying (Tt a )- Quite naturally, we construct a family of approximating nondegenerate problems 
as follows: Setting I2 r := 12 \ B r ( 0), we consider the boundary value problems 


(P)r 


f — div(cr(a;)Vu) = A u — \u\ 21 u, in fl r , 
l u\dn r = 0. 


From the characterization ifOll . problems (P) r are non-degenerate, and it can be shown that for 
fixed r > 0, there exists a global branch of positive solutions (see Definition 12.(ill , by using Harnack 
type inequalities. The next step is to prove that the limit of the approximating family ( P) r , as 
r —y 0, preserves the same property, thus Theorem o 

When is unbounded Cl, Lemma 2.2 and Remark 2.3] implies that under (H^P), in addition 
to 11. dl) . there exists R > 0, such that B r (zi) C B^(0) for every i,... ,k and 
(1.6) (in) a(x) > for x G fl \x\> R. 


In the unbounded domain case we consider as a model, the diffusion coefficient 
(1.7) a(x) = \x\ a + \x\ f \ a G (0, 2), /3 > 2, 


satisfying (JPp). Note that since a is unbounded the Harnack inequality is still not applicable. 
To approximate o defined in the unbounded domain (II C M. N ), this time we consider the 
approximate family of boundary value problems in flu := Q D Br( 0): 


(P)r 


f —div(cr(x)Vu) = Au — jul 2j u, in 

l = 0. 


Theorem o holds for (P)r and the claim is that as R —> oo the theorem remains valid at the 
limit. 

To establish the properties of the principal eigenvalues corresponding to both of the approximat¬ 
ing problems (P) r and ( P)r , we prefer an alternative proof, based on an appropriate adaptation 
of Picone’s Identity. This identity has been used in where the author established certain 

properties of the principal eigenvalue of the p-Laplacian operator, and extends Sturm Theorems 
to degenerate elliptic equations. 

Furthermore, we note that the presented method is applicable independently of the shape of H. 
In general, the situation becomes more complicated for non-uniformly elliptic problems in terms 
of u. As an example of the appearance of local bifurcation , for such a type of equation, we refer to 

m- 
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A general treatment of degenerate elliptic equations is provided by the monograph [2f)[ , focusing 
on the existence and properties of solutions (the issue of global bifurcation in the degenerate case 
is not addressed). Especially in the unbounded domain case, the problems are non-degenerate 
(at least in the sense of degeneracy, imposed by assumption (II.Old . In m a global bifurcation 
result is proved for a degenerate semilinear elliptic equation, with a degenerate diffusion coefficient 
of ’’critical exponent” (inducing non-compactness). Recent global bifurcation results for non¬ 
degenerate problems are included in the works 0 E2 EHUm EB1 ESI EZj For an overview, we also 
refer to the latest monographs fT5t . O). 

It is our intention to use Theorem o as a main tool, for a more detailed description of the 
asymptotic behavior of solutions of o, at least for the case of a bounded domain. A consequence 
of Theorem o is that for any fixed A > Ai, the set £ includes the trivial, the unique nonnegative 
solution of o and its (unique) nonpositive reflection. A combination of Theorems II. 111.51 could 
be used to design an intuitive picture for the dynamics of Ob It seems that the system undergoes 
through Ai a pitchfork bifurcation of supercritical type, where exchange of stability holds, i.e., the 
trivial solution is stable when A < Ai, while for Ai < A the nonnegative (nonpositive) solution of 
the global branch become the stable stationary state. Section 5 is devoted to some remarks related 
to the rigorous verification of the bifurcation picture for O- The fact that solutions of o with 
nonnegative initial data, remain nonnegative for all times (a ’’maximum principle” property), and 
the stability analysis of the unique nonnegative steady state, in conjuction with fTTH Theorem 2.7], 
implies the following 

Corollary 1.3. Assume that condition (H a ), holds. If (f>o > 0 a.e in f l, any solution <p(t) of ili.il) . 
tends to either the trivial or the unique nonnegative equilibrium point, as t —> oo. 

As it is expected, the nonnegative steady state is a global minimizer for the Lyapunov functional 
rRema.rk l5.3ll . A comment on the role of the “degeneracy exponent” a and a discussion concerning 
some possible further developments with respect to the case of noncompactness, is given in Remarks 

15.4115.51 

We conclude by mentioning the main results, on the convergence of globally defined and bounded 
solutions of evolution equations to rest points, as t —* oo. 

For scalar parabolic equations we refer to [551 611 and m for convergence to a single equilibrium. 
In (57) the result is proved for a semilinear heat equation defined in a higher dimensional domain, 
assuming a special structure of the set of rest points (semistable solutions). In convergence 
to a unique rest point, at least for the scalar case, is proved without the hypothesis that the set of 
rest points is totally disconnected. The same result is extended to semilinear parabolic and wave 
equations considered in multidimensional domains in 1-1511~T7T1 I I81 . when the nonlinearity is analytic. 
For a scalar degenerate parabolic problem (porous medium equation) a positive answer is given in 
m In the recent work the result of convergence to a (single) equilibrium is extended to a 
semilinear parabolic equation in 1 N : The main difficulty in the unbounded domain case is that 
even there exists a unique rest point z (radial with respect to 0), the w-limit set, may contains 
infinite many distinct translates of z. The authors introduce a new method, by defining moments 
of energy, which can discriminate against different translates of a rest point. The work m provides 
also a brief but complete review of the existing results and methods. For a more detailed survey 
we refer to Ell¬ 
in the case of non-autonomous systems or in the case where uniqueness of solutions of the 
evolution equation is not expected, the question on the convergence of solutions to rest points, and 
generally, on the existence of a global attractor, is discussed through the framework of generalized 
processes and semiflows in nnnana. Applications include nonautonomous semilinear wave 
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and parabolic equations, or equations involving non-Lipschitz nonlinearities, and Navier-Stokes 
equations. 


2. Preliminaries 


Function spaces and formulation of the problem. We recall some of the basic results on 
functional spaces defined in m■ Let -/V > 2, a G (0,2) and 

£e(2,+00), if a G (0,2), N = 2, 


2 =• 

"/-v • 


2 N 


( 2 ,^ 5 ), if a €(0,2), N >3. 


V-2+0! 

The exponent 2*, has the role of the critical exponent in the classical Sobolev embeddings. The 
following Caffarelli-Kohn-Nirenberg inequality holds, for a constant c depending only on (3, N, 


( 2 . 1 ) 


dx 


'r n 


< c 


\x\^\\7(j)\ 2 dx, for every 1 f G ^“(l^). 


By using CH> and conditions (H a ) and (H)f), it is proved in m Proposition 2.5], the following 
generalized version of cn>, 


( 2 . 2 ) 


\(j)\ 2a dx 


< K 


a(x)\V<f>\ 2 dx , for every cj) G Co°(f2). 


As a consequence of and eh , we have the following generalized version of Poincare 
inequality (ED Corollary 2.6-Proposition 3.5], see also J3J3 Section 5]). 


Lemma 2.1. Let LI be a bounded (unbounded) domain of 1^, N > 2 and assume that condition 
(Ha) ((H)))) is satisfied. Then there exists a constant c > 0, such that 

(2.3) f |<^>| 2 cZa; < c f er(x)|V<)>| 2 d:r, for every </> G CjJ°(f2). 

J Q, JQ 

We emphasize that inequalities EH , EH (and EH in the case of a bounded domain), hold 
for some a G (0, 2], However, the case a = 2 can be considered as a “critical case” with respect to 
compactness of various embeddings , even in the bounded domain case. Moreover, condition (H a ) 
is optimal in the following sense: For a > 2 there exist functions such that EH is not satisfied 
m- Note also that in the case of an unbounded domain, EH does not hold in general, if (3 < 2 
in (H)f). We refer also to the examples of p. 

The natural energy space for the problems EH and EH involves the space Uq' 2 (H,ct), defined 
as the closure of Cq°(LI) with respect to the norm 

■= (^(*) l^l 2 ) • 

The space D)’ 2 (fl, a) is a Hilbert space with respect to the scalar product 

(</>, ip) a := / a(x)V^Vt/j dx, for every <j> if G Dl’ 2 (Ll, <j). 

J n 

The following two lemmas refer to the continuous and compact inclusions of D)’ 2 (LI, a) ) T5S 
Propositions 3.3-3.5]. 

Lemma 2.2. Assume that LI is a bounded domain in 1*, N > 2 and a satisfies (H a ). Then the 
following embeddings hold: 

i) D)' 2 (Ll, a) L 2 <* (12) continuously, 

ii) D)' 2 (Ll,a) L p (Ll ) compactly if p G [1,2*). 
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Lemma 2.3. Assume that Cl is an unbounded domain in R^, N >2, and a satisfies (Tt^f). Then 
the following embeddings hold: 

i) Dq' 2 (CI, a) L p (Cl) continuously for every p G [2£, 2*], 

ii) Dq’ 2 (CI,ct) L p (Cl) compactly if p G (2^,2*). 

Remark 2.4. It is crucial to note that as a special case, the embedding Dq’ 2 (Q, a) C L 2 (Cl) is 
compact if either conditions (TL a ) or (Tt^) hold: Observe that (3 > 2 implies 2 ^ = N ™ + g < 2 , 
i.e. 2 G (2^,2*). In the unbounded domain case, we need a to grow faster than quadratically at 
infinity, to ensure compactness. We also stress the fact, that since a is not in L^ c (Cl), there is 
not in general any inclusion relation between the space Df (Cl, a) and the standard Sobolev space 

H l(Ci). 


To justify the natural energy space for o, we have applied in ESI, Friedrich’s extension 
theory ED pg. 28, 32], [73 pg. 126-135]: Assuming conditions (TL a ) or (7fg°), the operator 
T = —div(cr(a:)V(/>) is positive and self adjoint, with domain of definition 

D(T) = [ct> G Dl'\Cl, a), T0 G L 2 (H)} . 

The space D(T), is a Hilbert space endowed with the usual graph scalar product. Moreover, there 
exist a complete system of eigensolutions {ej, Ay }, 


z 24 ) f -div(cr(a:)Vej) = Aj-ey, J = 1,2,..., j> 1, 

' \ 0 < Ai < A 2 < A j —> 00 , as j —► 00 . 

The fractional powers are defined as follows: For every s > 0, T s is an unbounded selfadjoint 
operator in L 2 (Cl), with domain D(T S ) to be a dense subset in L 2 (Cl). The operator T s is strictly 
positive and injective. Also, D(T S ) endowed with the scalar product (cj),if) d(t s ) = (T s 0, T s -0) i 2, 
becomes a Hilbert space. We write as usual, V 2 s = D(T S ) and we have the following identifica¬ 
tions D(T -1 / 2 ) = D^(CL,a)= the dual of D^ 2 (Cl,a). D(T°) = F 2 (Cl) and D(T 1 /2) = D^ 2 (Cl,a). 
Moreover, the injection V 2 Sl C V 2 s 2 , s i, s 2 GR, si > S 2 , is compact and dense. 

While in EH, the local in time solvability was discussed via compactness methods, for the 
purposes of the present work, it is more convenient to study the local in time solvability of o 
in Dq (Cl, a), via the semigroup method approach: The discussion above clearly shows, that the 
operator — T is the generator of a linear strongly continuous semigroup T(t) ( | l()l 1221163| h 


Definition 2.5. For a given function <fo G Dq’ 2 (CI, a), 0 < 7 < 00 and T > 0, a solution for the 
problem D. is a function 

<t>(x, t ) € C([0, T); Dl' 2 (Cl, a)) n C^O, T); L 2 (H)), 


satisfying the variation of constants formula 


(2.5) 


<t>(t) = T(t)<j> 0 + / T(t - s)f(<j>(s))ds 


where f(s) = As — |s| 27 s. 


Solutions of HU satisfying Definition 12.51 and solutions satisfying 2T>. Definition 2.3] (weak 
solutions) are the same. This is an immediate consequence of 

We conclude this introductory section, by stating for the convenience of the reader, some basic 
definitions and results for our analysis. We state first a result on the existence of a branch of 
solutions of an operator equation (bifurcation in the sense of Rabinowitz m -see also ESI). 
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Theorem 2.6. Let X be a Banach space with norm || • ||a' and consider the operators 

G(A, ■), L, H(A, ■) : X —¥ X*, 

where G(A, •) = AL(-) +H(A, •), L is a compact linear operator and H(A, •) is compact and satisfies 

lim PMk =l 

IMIx-o \\u\\x 

If X is a simple eigenvalue of L then the closure of the set 

C = {(A, u) G R x X : (A, it) solves N(A, u) := u — G(A, u) = 0 in X*, u ^ 0}, 

possesses a maximal continuum (i.e. connected branch) of solutions, C\, such that (A, 0) G C\ and 
C\ either: 

(i) meets infinity in Mxl or, 

(ii) meets (A*,0), where A* ^ A is also an eigenvalue ofL. 

In the approximation procedure, we are making use of a generalized Harnack-type inequality 
(see |2H1 HT7] and the references therein). 

Theorem 2.7. (Harnack-type Inequality) Consider the equation 
(2.6) — div(a(x, u)\Vu\ p ~ 2 Xu) = f(x, u), x G Cl, 

where f l C M. N , 1 < p < N and the functions a and f satisfy the following conditions: 

(i) a is a Caratheodory function, such that a(x,s ) is uniformly separated from zero and bounded 
for almost every x G LI and all s G M, 

(ii) f is a Caratheodory function and for any M > 0 there exists a constant cm > 0, such that 

|/(x,s)| < c m \s\ pt -\ 

for almost every x G Cl and all s G (— M,M), where p* is the critical Sobolev exponent p* = . 

Assume that u G £> 1 ,p (n) := {u G L p (f2) : Vm G (L p (fl)) N } is a weak solution of 12.6 1) satis¬ 
fying the weak formula 


a(x, it) | Vrt| p 2 XuS7(j)dx= / f{x,u)4>dx, 
! Jn 


holds for any (j> G Co°(f2). Then, for any cube K = K(3p) C ft with 0 < u < M in K, we have 
that 

max u(x) < C min u(x). 

x£K p x&K p 

In particular, if the weak solution u ^ 0 of m satisfies u > 0 in then it follows that u is 
strictly positive in Cl. 

Remark 2.8. In the case where a{x,s) = a{x) satisfies condition (i) of Theorem 12. 61 the norms 
of Dl' 2 (Cl, a) and D 1 ’ 2 (Cl) are equivalent. 

We also recall some basic definitions and results on semiflows (see U2II3] and Ell ED)- Let x 
be a complete metric space. For each </> 0 G X, via the correspondence S(t)<j > 0 = <j>(t), a semiflow is 
a family of continuous maps S(t) : X — > X, t > 0, satisfying the semigroup identities (a) 5(0) = I, 
(b) S(s + t) = S(s)S(t). For B C X, and t > 0 

S(t)B := {4>(t) = S(t)(f >o with ^>(0) = <po G B}. 

The positive orbit of <j) through </> 0 is the set 7 + (^o) = {<t>(t) = 5(t)</>o, t > 0}. If B C X then the 
positive orbit of B is the set 

7 + (B) = (J S(t)B = { 7 + ( 0 ) : <j>(t) = 5(f)0 o with ^( 0 ) = (fo G B}. 

t> 0 



If to > 0, 7 ta (B) := U t>t0 S(t)B = j + (S(to)B). The w-limit set of fio £ X is the set w(</>o) = {z G 
X : = S(tj)4>o —■• 2 for some sequence tj —* +oo}. A complete orbit containing <p 0 £ X , is 

a function </> : R —> X such that </>(0) = <po and for any s € R, S(t)(f>(s) = <j>(t + s ) for t > 0. If ^ is 
a complete orbit containing <p 0 , then the a-limit set of <j )o is the set 

a(<j> o) = {z £ X : 4>{tj) —> z for some sequence tj —> — oo}. 

The subset A attracts a set B if dist(«S(i)B, A) —> 0 as t —> +oo. The set A is positively invariant 
if S(t)A C A , for allt > 0 and invariant if S(t)A = A for all t > 0 . The set A is a global attractor 
if it is compact, invariant, and attracts all bounded sets. 

The semiflow S(t) is eventually bounded if given any bounded set B C X, there exists f 0 > 0 
such that the set 7 to {B) is bounded. The semiflow S{t) is said to be point dissipative if there is 
a bounded set Bq that attracts each point of X. It is called asymptotically compact if for any 
bounded sequence 4> n in X and for any sequence t n —► 00, the sequence S(t n )(f> n has a convergent 
subsequence. It is called asymptotically smooth if whenever B is nonempty, bounded and positively 
invariant, there exists a compact set K. which attracts B. 

A complete orbit is stationary if cf>(t) = z for all t £ R for some z £ X and each such z, is called 
an equilibrium point. We denote by £ the set of stationary points. 

The functional J : X —> R is a Lyapunov functional for the semiflow S(t) if (i) J is continuous, 
(ii) J{S{t)(j)o) < JS{s)(j) 0 ) and t > s > 0, (iii) if J{<j>(t)) =constant for some complete orbit </> 
and all t £ R, then </> is stationary. 

To derive the convergence result we shall use the following Theorem. 

Theorem 2.8. /|lli|) Let S(t) be an asymptotically compact semiflow and suppose that there 
exists a Lyapunov function J. Suppose further that the set £ is bounded. Then S(t) is point 
dissipative, so that there exists a global attractor A. For each complete orbit containing <fio lying 
in A the limit sets a((f>o) and u>(cj) 0 ) are connected subsets of £ on which J is constant. If £ is 
totally disconnected (in particular if £ is countable) the limits 

z- = lim z+ = lim <f>(t) 

t —► — OO t —» + oo 

exist and are equilibrium points. Furthermore any solution S(t)<f 0 tends to an equilibrium point as 
t —> 00 


3. Global Attractor in D l f 2 (£t,a) 

In this section we shall show, that the degenerate semilinear parabolic equation o> defines a 
semiflow in the energy space D(f (f l, cr), possessing a global attractor. We state first an auxiliary 
lemma. 


Lemma 3.1. Assume that either conditions (Tt a ) or (Tip 3 ) hold. The function f\(s) 
R, defines a sequentially weakly continuous map fi : T>Q(£t,cr) —> L 2 (fl) if 


(3.1) 


0 < 7 < 


2-a 


2(N-2 + a) 


:= 7 


|s| 27 s, s £ 


Furthermore, if Fi(<p) := fi(s)ds, the functional Ei : Dq’ 2 (Q,(t) —> R defined by = 

f n Fifjfjdx, is C' 1 (Z?Q , 2 (n, <t),R) and sequentially weakly continuous. 


Proof: It can be easily checked that the functional fi is well defined, under the restriction 
EH), by using Lemmas 12 .2f il- l!OT ih Similarly, it follows that Ei is well defined if 

2-a 


0 < 7 < 


N-2 + a 
9 


■■= 7i) 


(3.2) 




and note that 7 * < 71 . To show that both functionals are sequentially weakly continuous, we use 
the compactness of the embeddings stated in Lemmas m iil- l2.3f ii). and repeat the lines of the 
proof of El Lemma 3.3, pg. 38 & Theorem 3.6, pg. 40]. To verify that Ei is a C 1 -functional, and 
its derivative is given by 

(3.3) E[((t>)(z) = (fi(4>),z) , for every 4> £ Dq’ 2 (CI,<j), 2 : e D(( 1 (f], a), 


we consider for £ ZDq’ 2 (O, ct), the quantity 

Ei(4> + sip) - E i(<^>) 
s 

(3.4) 


-If - 777^1 (0 + 9sip)d9dx 
s Jn Jo d0 

f I fi(cf> + s9il’)d6dx. 

J o Jo 


Setting q = N f_ 2 _ a ’ 9 + 2* = 1, we observe that 


(3.5) 



+ 9sip)ipdx 



(27+1) 


\ip\ ( 27 + 1 ) ) 9 d^ 



Lemmas 0 i) 0 i) are applicable under the requirement (2y + 1 )q < 2* which justifies (13.21) . 
Using the dominated convergence theorem, we may let s —> 0, to obtain that E is differentiable 
with the derivative E!• 

We consider next a sequence {(j4 n }neN of D 0 ’ (U, a) such that <j> n —* <f> in D 0 ’ (f l, a) as n —> 00. 
It holds that 


(3.6) 


(E'M n ) - E[(<t>), z) < || h{<t>n) - /l(^)||£,|klfe 


2 * 

Setting pi = ^ we observe that the requirement for p 1 > 1, justifies the restrictions on the 
exponent of degeneracy a , imposed by (Tt a ) or ("^p 3 )- Setting now p 2 = 0° 1 P 2 1 + Pf 1 = 1, 

we get 

ll/i(<An) - fi{4>)\\li < + \<t>\ 2l ) qP2 dx^j (X 1^™ ~~ ^I 2 “* dx ) 


Let P 3 = 2^qp 2 . To apply Lemmas E3ti)Di) once again, we need P 3 < 2 * or l13.2[) . 

Under this condition we have that A (4> n , 4>) —> 0 as n —* 00 and from 113.61) . we get the continuity 
of E'. o 

We consider the energy functional J : Z?q’ 2 (U,ct) —> M 

(3.7) J ^ cr(x)\S/(/>\ 2 dx - ^ \(p\ 2 dx + ^ J ^ \(p\ 2 l+ 2 dx. 


Proposition 3.2. Let <j> 0 £ Dq’ 2 (U, a) and either conditions (Tta) or (TEp) be fulfilled, and assume 
that Q holds. Then equation D. has a unique, global in time (weak) solution (f), such that 

( 3 . 8 ) cf £ C([ 0 ,00); Dl’ 2 (n, a)) n C^O, 00); L 2 (U)). 

For each (weak solution) (</>(■)) £ C 1 ( [ 0 , 00)) and 

( 3 . 9 ) ^ J (</>(f )) = ~ J \dt<f>\ 2 dx 


10 








Proof: By using similar arguments to those used for the proof of Lemma I.'!. 1 1 we may show 
under the assumption El, that the function /(s) = |s| 27 s — As, defines a locally Lipschitz map 
/ : Dq (fi; a) —> This suffices in order to show the existence of a unique solution <fi with 

</>(0) = cj)Q , defined on a maximal interval [0,T mox ), where 0 < T max < oo [22]. 

We proceed by showing that T max = oo. First note, that by Lemma ld.il the energy functional 
J is C 1 . This fact allows to adapt the method of PH mu, in order to justify Id.91) for any t G [0, T], 
T < T max . We repeat the main lines of the proof, only for the shake of completeness: For all 
<fi € D(T), then 


<-T0 + f (</>), = - I \&w{a{x)V4>)+fm 2 dx 

Jn 

(3.10) = — f \d t (j)\ 2 dx < 0. 

Jn 

Setting g(t) = f((j){t)) we consider sequences g n {t) G C 1 ([0, T]; Dq' 2 (Q.\ a)) and </>o n G D(T) such 
that 


9n -> g , in C 1 ([0,T]-Dl' 2 (fl-a)), 

4>0n -*■ 4>o, in Do' 2 (fl;cr). 

We define <j> n {t) = T(t)(t> 0 n + /q T{t — s)g n (s)ds 1 and it follows from [B3J Corrolary 2.5, pl07] that 
4>n{t) G D(T), cj) n G C 1 ([0, T]; Dq’ 2 (Q; a)) satisfying ^</>„(f) + Tcj> + /(</>) = 0. Also, from ;Il :l 
Lemma 5.5, pg. 246-247] (see also [I3J Theorem 3.6, pg. 41]) we get that 

in C'([0,T];Do’ 2 (fl;cr)). 

Now using the fact that J is C 1 and font . we may pass to the limit to 

~ J{4>0n) = [ (i7'(<Ms)), -T0„(s) +g n (s))ds 

Jo 

= - f \\dt ( f > n{s)\\i J 2 ds + [ {J'{4> n {s)),g n {s) - f((j> n {s)))ds 

Jo Jo 

to derive EH]). 

Multiplying EB by (j), and integrating over f 2 , we obtain the equation 


(3.11) ~||*||i.+ [ (j{x)\V<t>\ 2 dx- XM \ 2 L 2 + f \4>\ 2l,+2 dx = 0. 

We are focusing on the case where A > Ai and the domain is unbounded. By interpolation and 
Lemma IP fib we have that for some 6 G (0,1), 




< 


-e) 


< 

2XCf~ e) \ |^|||V2 

m 

,2(1-0) 

< 

1 

2 

ll^llDj' 2 (f2,<T) +Cl 

m 

L 2 7+ 2 

< 

1 

2 


m 

27+2 , 77 

£27 + 2 T JT-l 
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(3.12) 


where Cp is the constant of the embedding Z}°' 2 (fi, a) C L 2 p(£1) and 

(7+1)08-2) „ • 


0 = 


(AT + /? — 2)(q + 1) — AT ’ 


ci = C$(2 A)t=* (1-0) (20)i 


(3.13) 


i?l = C : 


2±l 2-V7 


(7 + 1 )T +1 ’ 

By inserting the estimate cmt to ehj, we get 

2^ll0lli 2 + 2H^Hz)J' 2 (n,<T) + + 911011^27+2dz < R-i- 

Gronwall’s Lemma leads to the following inequality 

(3.14) \\4>(t)\\ 2 L* < ||^(0)||l a exp(—2At) + —!^-(l - exp(-2At)). 

Letting t —► 00 , from etta we obtain that 

(3.15) limsup || 0 (t )|| 2 2 < p 2 , p 2 = 1/ARi. 

t—* OO 

Now assume that <j>o is in a bounded set $ of Dg’ 2 (fl, cr). Then (13.151) implies that for any p\ > p , 
there exists t 0 (B,p ±), such that 


(3.16) || 0(£)]| L 2 < p 1; for any t > to(B,pi). 

We observe that by the definition of the energy functional and (EHTH) . 

J{<!>{*)) > \ f <x{x)\S7 (j )\ 2 dx — ^ f \(j>\ 2 dx 

(3.17) > - j^a(x)\V(j )\ 2 dx - ^p\, t>t 0 . 

Hence, since £f((f>(t)) is nonincreasing in t, we conclude that 

(3.18) H^( i )lli3^ 2 (n,<j) — 277(0°) + \p\, t>t 0 . 

Thus solutions are globally defined in D°’ 2 (fi, cr). o 

Proof of Theorem o It is not a loss of generality to assume that 4>o G B(0,R), a closed 
ball of Dq (0, cr), of center 0 and radius R. Then from Lemma 13. II and (13.71) it follows that there 
exists a constant c(R) such that J(4>o) < c(R). Hence, 13.1811 implies that S(t) is eventually 
bounded. Since the resolvent of the operator — T is compact, S(t) is completely continuous for 
t > 0, thus asymptotically smooth. The equivalence criterion m Proposition 2.3, pg. 36], implies 
that S(t) is asymptotically compact. The positive orbit 7 + (</>o) is precompact, having a nonempty 
compact connected invariant w-limit set u>((j> 0 ). From (13.611 and the continuity of S(t) it follows 
that G £. 

It remains to show that £ is bounded, to conclude that S(t) is point dissipative. An equilibrium 
point of S(t), is an extreme value of the functional J or equivalently, satisfies the weak formula 


(3.19) / a{x)VvS7vdx — A / uvdx + / |it| 27 icu<&E, for every v G Z?°’ 2 (fl, cr). 

if! is! Jfl 

Setting v = u in (EH31) and using inequality E321) we obtain 


(3.20) 


<r(x)\S7u\ 2 dx < (2A) to-®) i? 2 ( 7 , 0), 


(see nm which implies that for fixed A the set £ is bounded, o 
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4. Global Bifurcation of Stationary States 


The validity of the continuous imbedding L>q’ 2 (D, a) L 2 <*(fl) ('Lemmas 12.212.31 fill, enables 
us to use the same arguments as in the proof [2TH Theorem 4.1 (Step 2)] (see also P Lemma 2.8]), 
in order to prove L°°-estimates, for the weak solutions of o and JH. 


Lemma 4.1. Assume that Q is an arbitrary domain (bounded or unbounded) and the conditions 
( , H a ) or (7f“) are satisfied. Then any weak solution u of ill.HI) or El is uniformly bounded in 
fl, i.e. ||u||j,oo(q) < C, where C is a positive constant depending on A, 7 and K, where K, is the 
constant appearing in El - 


A. The bounded domain case: We assume that the diffusion coefficient is given by d, 
and we consider the following problems: 


(P) 


( — div(u(a;)VM) = A u — \u\ 21 u, in SI, 

l u\on = 0, 


( — div(tr(x)Vu) = Xu — |u | 27 u, in f l r = fl\B r ( 0), 

w) 1 " 1 1 n 

l = 0, 

for some r > 0 sufficiently small. Standard regularity results (cf. pTTTl Theorem 8.22]) imply that 
if u is a weak solution of the problem (P), ((P) r ) then u £ G i 1 o ’^(fl\{ 0 }), (u £ C^(Cl r )), for some 
( 6 ( 0 , 1 )- 

For the linear eigenvalue problems 

—div(<r(a:)Vu) = A u, in Q (fl r ), 

w|an = 0, (M|an r = 0), 

we have the following lemma. 


(PL), ((PL) r 


Lemma 4.2. Assume that a is given by El- Problem (PL) ((PL) r ), admits a positive principal 
eigenvalue Ai (Xi jr ), given by 


(4.1) 


Ai (Ai, r ) — 


inf 

<j> £ Dq' 2 (XI (fl r ), cr) 


In (n r ) a ( x ) l V( ^l 2 dx 
In (n r ) 1^1 2 dx 


with the following properties: (i) Ai (Ai ir ), is simple with a nonnegative (positive) associated 
eigenfunction u\, (ui r ). (ii) Ai (Ai ir ), is the only eigenvalue of (PL) ((PL) r ), with nonnegative 
(positive) associated eigenfunction. 

Proof: The existence of Ai (Ai, r ) is a consequence of Lemma [2.21 (fif (see also (12.41) 1. For the 
proof of (i), let us assume that iti > 0 (iti jr > 0) in S l (f2 r ) (since if u (wi,r) is a minimizer of (14.11) . 
then |mi| (|ui ir |) must be also a minimizer-similar arguments may also be find in |39l Theorem 
8.38]). The simplicity of Ai (Ai iT .), can be shown by an alternative argument, based on the so 
called Picone’s Identity PI3I1 

(PX): Assume that u > 0, v > 0 are almost everywhere differentiable functions in SL Define 


L(u, v) 
R(u , v) 


2 

|Vu | 2 + ^|Vu| 2 -2-Vu-Vw, 

v z V 

|Vu| 2 - ■ Vu. 
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Then L(u,v) = R(u,v), L(u,v) > 0, and L(u,v) = 0, if and only if u — kv for some constant k, 
a.e. in d. 

Let flo C Cl a compact subset of d, and 0 < (f> £ Cg°(d). For A > 0 we consider u G C]f c (d), 
( 6 (0,1), a weak solution of ( PL ), such that 0 < u a.e in Cl. Then, for any e > 0, we have that 


0 


(4.2) 


/ a(x) L(cf>,u + e ) dx < / a(x) L((f,u + e) dx = 

J Qq J f2 

/ a(x) R(4>, u + e) dx = 

Jn 

[ a(x) \S7(j)\ 2 dx — [ a(x)V\ — — J • Vm dx = 

J n J n \ u + e J 

f a(x) \S7(j)\ 2 dx + [ (— —)V(<7 (x)Vu)dx = 

Jn Jn\u + eJ 


f a{x) |V0 | 2 dx — A f 
Jn Jn 


u + e 


u dx. 


Assume now that Ai is not simple. Let v ^ U\ be another associated eigenfunction, v G 
Z?Q , 2 (n,cr) almost everywhere differentiable in d, such that v(x) > 0 in some fl + C d. Consider 
iTOll with do C A = Ai and u = u\. Letting cf> v in and e — * 0, Fatou’s Lemma and 
Lebesgue Dominated Convergence Theorem, imply that L(v,u\) = 0 a.e. in f2 + . Hence from (VI) 
we get that v = kui, a.e. in d + , which implies the simplicity of Ai. Property (i) is proved. 

For the proof of (ii), we suppose that there exists another eigenvalue of (PL), A* > Ai, to which 
corresponds a nonnegative eigenfunction u*. Consider (TOl) with do C A — A and u — u*. 
Letting <j> —* u\ in d and e —> 0, we obtain that 


which is a contradiction, o 


0 < / a(x) L(u\,u*) dx < 0, 

Jn 


Lemma 4.3. Assume that a is given by Let also Ai, Ai >r , be the positive principal eigen¬ 

values of the problems (PL), (PL) r , respectively. Then, uq r —» u\ in D^’ 2 (d, a) n L™ c (d \ {0}), 
and Ai jT . | Ai, as r J. 0. 


Proof: We extend ui, r on d as 


ui,r(x) =■■ 


Ui, r (x), X G d r , 

0, x G B r , 


for any sufhciently small r > 0 , but in the sequel, for convenience, we shall use the same notation 
u\ t r = Ui, r . Observe that 


Ai r — 


fn r a ( x ) |Vui >r | 2 dx J n cr(x) |V-ui , r | 2 dx 


In, Krl 2 dx 


fn \ u Lr \ 2 dx 


> Ai, 


and Ai jT - is an decreasing sequence, as r —> 0, since d p C d g , for any p > g . Clearly, u i, r forms 
a bounded sequence in Dq’ 2 (H, a). Tjemma f 2 . 2 l (ii) and Lemma roi imply the existence of a pair 
(A *,u*), and a subsequence of u\, r (not relabelled), such that 


Ai- 


\u\ r Ydx —* A* / \u*\ 2 dx 
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as r —► 0. Then <EU) implies that u\^ r —> it* in Z)q’ 2 (0, ct) and u* satisfies 

f <y(x) |Vw *| 2 dx = A* f \u*\ 2 dx. 

J n J n 

From Lemma E 21 (ii), we obtain that (A *,u*) = (Ai,iti). We conclude by justifying the claim 
that U\, r is uniformly bounded in Note that Ayr £ (Ai, Ai + e), for some e > 0 and any r 

small enough. Since ui, r £ -Dq’ 2 (S2, ex), it holds that ||wi,r||iaj(n) < ^ W^xWd 1 ’ 2 ^ ^ is gi yen 

in (l 2 ~ 2 l) and is independent of r. Hence, from Lemma rra we have that iti. r is uniformly bounded 
in L°°(fl). Then by a standard bootstrap argument we get that iti, r —> u\ in L^ c (fi \ {0}) and the 
proof is completed, o 


Proposition 4.4. Assume that a is given by Ob The principal eigenvalues \±, Ay r of the linear 
problems {PL), {PL) r , are bifurcation points of the problems (P), ( P) r respectively. Moreover, 
for any (sufficiently small) r > 0, the branch C\ l r is global, and any function which belongs to 
C \i , is strictly positive. 


Proof: The existence of branches bifurcating from Ai, Ai, r follows by Theorem !!!, (il since Lemma 
12.21 (ii) and Lemmas 14.214.HI are in hand. We outline the proof for the branch C\ t . 

As in W, we define a bilinear form in C , Q°(fl) by 

(4.3) (u, v) = [ cr{x)S7uVv dx — f uvdx, for all u,v £ C^°(f2). 

J n 2 Jq 

(c is the constant in 12.31 1 and we define A' to be the completion of Cq°(Q) with respect to the 
norm induced by 14. 3D . ||u[|^- = (u,u): from inequality 12.31 we get that 


1 

2 


'PT(SV) 


< 


Mix 


< 


'Df 2 (n,ay 


for all u,v £ 


and by density it follows that X = Z?q’ 2 (H, ct). Henceforth we may suppose that the norm in X 
coincides with the norm in Pq’ 2 (H, a) and that the inner product in X is given by < u,v >= 
( u,v) a (moreover, we may assume that if < •,• >.y,x* denotes the duality pairing on X , then 
< ’) ’ >X,X*=< ' > £ 3 - Identification Principle 21.18, pg. 254]). On the other hand, the bilinear 

form 


a (it, v) = / uv dx, for all u, v £ X, 

Jn 

is clearly continuous in X as it follows from Ijemma. l 2 . 2 l . and by the Riesz reperesentation theorem 
we can define a bounded linear operator L such that 


(4.4) &{u,v) =< L u,v >, for all u,v £ X. 

The operator L is self adjoint and by Lemma o (ii) is compact. Th largest eigenvalue v\ of L is 
given by 

< L u, u > Jo u 2 dx 

Vl = sup —-— = sup - - - . 

utx <u,u> ueX J n a{x)\\7ufdx 

It follows from Lemma FOl that the the positive eigenfunction ui of {PL) corresponding to Ai is a 
positive eigenfunction of L corresponding to = 1/Ai. We consider now the nonlinear operator 
N(A, •) : R x X —> A'* defined by 


(4.5) (N(A, it), v) = / <t{x)XuXv dx — A / uvdx+ / \u\ 2l uvdx, for all v £ X. 
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Arguments very similar to those used for the proof of Lemma Id. II can be used in order to verify 
that for fixed u £ X, the functional S defined by 

S(v)= [ cr(x)XuVvdx- X f uvdx+ [ \u\ 2 l uvdx, v £ X, 

is a bounded linear functional and thus N(A, it) is well defined from (14.511 . Moreover by using 
the fact that X = D ( ]' (S~2, cr) and relation 14.411 . we can rewrite N(A, u) in the form N(A, it) = 
u — G(A, u) where G(A, u) := ALit — H(it), 

<H(it),i;>= / \u\ 2l uvdx for all v £ X. 

J Q 

The restriction ED and Lemma E2 (ii) implies that H is compact. Moreover we observe that 


IMIx 

(4.6) 

Therefore, we get from (IQ1) that 

'|H(u)||x* 


< H(it), v > | < 


\u\\ x 


1^27+2 I M |l 2 2'+ 2 1 M | L 2 "v+ 2 


< Cl lilt 




V \x- 


lim 

Mix—>0 


= lim 


1 


sup 


< H(it), v > | = 0 . 


.Mix IMIx-»o|| 0 || x <i ||w||x 
To prove that C\ l r is global for sufficiently small r > 0, we proceed in two steps. 

(a) We shall prove first that for all solutions (A, it) £ C\ liT close to (Ai )r ,0) it holds that u(x) > 0, 
x £ f l r . In other words, we have to show that there exists eo > 0, such that for any (A,it(x)) £ 
C\ l r fl B eo (( A 1t , 0)), it holds that u(x) > 0, for any x £ Sl r (By B eo ((Ai ir , 0)), we denote the open 
ball of C\ l r of center (Ai, r ,0) and radius eo). 

We argue by contradiction: Let (X n ,u n ) be a sequence of solutions of (P) r , such that (X n ,u n ) —> 
(Ai /r , 0) and assume that u n are changing sign in fl r . Let u~ := min{0, u n } and U~ =: {x £ tt r : 
u n (x) < 0}. Since u n = — u~ is a solution of the problem (P) r it can be easily seen that u~, 

satisfies (in the weak sense) the equation 

(4.7) - div (a(x)Vu~) - X n u~ + \u n \ 2l u~ = 0, 

u n law = 0 . 

Then, multiplying ED with u n and integrating over we have that 


(4.8) 


)Un 


|Vit„ \ 2 dx - X n 


—12 


dx - 


121 . — 12 


dx = 0 . 


>Un 


)Un 


Since A„ is a bounded sequence, it follows from EEL Holder’s inequality and relation ED that 

ll u «Hr>J’ 2 (w-,<r) - Xn I \ u n? dx 


c I u: 


U r 


>Un 

— C^A n \ N ll u nll D i.2 

or, equivalently 

(4.9) M < \U~\, for all n, 

where the constant M is independent of n. We denote now by u n , = itn/||itn|| the normalization 
of u n . Then there exists a subsequence of u n (not relabelled) converging weakly in £>Q , 2 (fl r , a) to 
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some function uq. By following the lines of the proof of Lemma 14.HI it can be seen that uo = «i, r - 
Moreover, u n —> «i ir > 0 in L 2 (Q r ). Passing to a further subsequence if necessary, by Egorov’s 
Theorem, u n —> ui :T uniformly on Q r with the exception of a set of arbitrary small measure. This 
contradicts ifOll and we conclude the functions u n cannot change sign (for a similar argument, 

we refer to EH1I23IS3)- 

(b) Suppose now that for some solution (A ,u) £ C Xl , r , there exists a point xq £ fl r , such that 
u(xq) < 0. Using (a), the fact that the continuum C Xljr is connected (see Theorem 12.(ill and the 
C]£ (Q r )- regularity of solutions, we get that there exists (Ao,wo) G C Xl ^ r , such that uq(x) > 0 , for 
all x £ fexcept possibly some point xo £ Q r , such that u o(xo) = 0. Then Theorem 12.71 implies 
that Mo = 0 on Q r . Thus, we may construct a sequence {(A„,u n )} C C Xl , r , such that u n (x) > 0, 
for all n and x £ f l r , u n —> 0 in Df (f 1 r , a), and A„ —> Ao- However, this is true only for Ao = Ai- 
As a consequence, we have that C Xir cannot cross (A, 0) for some A ^ Ai, and every function 
which belongs to C Xl>r is strictly positive, o 

Theorem 4.5. Assume that a is given by Then, C\ x is a global branch of nonnegative 

solutions for the problem (P). 

Proof It suffices to prove that C\ l r —> C\ 1 , as r —> 0. The global character of C\ l r implies 
that for any fixed positive number R , and any r sufficiently small, the set C\ l r fl Pfl(Ai j7 ., 0) is not 
empty. By using the properties of Ai established in Lemma 14.21 and the compactness arguments of 
Lemma FOl we can show that 

lirn C \ 1 r n B r { Ai, r , 0) -> C Xl n B r { Ai, 0), for every R > 0, 

r—> 0 

which implies that C Xl r —> C Xl , as r —> 0. Alternatively, one may use Whyburn’s Theorem 

[3113113211211- 

Proof of Theorem [EH in the case of (H a )' One has to extend Theorem 14.51 in the case 
of a diffusion coefficient satisfying (H a ). Since the set of zeroes of a, Z a is finite, we may use 
<EJ and consider approximating problems similar to (P r ), defined this time in the domain Q r = 
fi \ (Ji B r (zi). The finiteness of Z a , allows to repeat the proofs of Lemmas 14 . 214.HI and Proposition 
IP1 without additional complications. 


B. The unbounded domain case We assume that the diffusion coefficient is given by o and 
we consider the following problem: 


(-P)oo 


—div(cr(a;)VM) = A u — \u\ 21 u, in fl, 
w|dn = 0. 


where H C R w , N > 2, is an unbounded domain containing the origin. The regularity results 
of [22 Theorem 8.22], imply once again that if u is a weak solution of the problem (P 0a ), then 
u e C«(O\{0}), for some f £ (0,1). This time, we consider the approximating problem, 

—div(tr(a;)VM) = A u — |u| 27 u, in Qr = Q fl Br{ 0), 

«|dn R = 0 . 



We consider the linear eigenvalue problems 


(■ PL )oo, ((PL) r ) 


( — div(cr(:r)Vu) = A u, in Q (Qr), 
l u|an = 0, (w|dn R = 0). 


A result similar to Lemma rr~2i holds. 
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Lemma 4.6. Assume that a is given by O- Problem ( PL ) OQ ((PL)r), admits a positive prin¬ 
cipal eigenvalue Ai (Xi, R ), given by 


(4.10) 


Ai = 


inf 

e Dl’^n {n R ),a) 


U 


si (r2j 


a(x) |V0| 2 dx 


h 


n (On) 


|0| 2 dx 


0^0 


with the following properties: (i) \± (Xi,r), is simple with a nonnegative (positive) associated 
eigenfunction u\, (u\,r). (ii) Ai (A^r), is the only eigenvalue of(PL) oc (( PL)r), with nonnegative 
(positive) associated eigenfunction. 


To prove a similar to Lemma roi result, we shall use the extension 

«i,r(z), x £ Hr, 


ui ,r(x) =: 


0 , 


x £ LI\LIr, 


and use for convenience the notation u\r = m^r. 


Lemma 4.7. Let Ai, (\i,r) be the positive principal eigenvalues of the problems (PPoo), (( PL)r ). 
Then, Ui t R —> u\, in ZAq’ 2 (L2, a) fl \ {0}) and Air | Ai, as R —» oo. 

We remark that for each R > 0, Theorem run is applicable for (P)r: There exists a global 
branch, C\ 1R , of nonnegative solutions, bifurcating from A^r. This suffices for a repetition of 
arguments similar to those used for the proof of Theorem run to show that C\ 1R —> C\ 1 , as 
P —» oo. 


Theorem 4.8. Assume that er is given by O- Then, \\ is a bifurcating point of the problem 
(Poo) and C\ x is a global branch of nonnegative solutions. 

Proof of Theorem 11.21 in the case of (' H One has to consider approximating problems 
similar to (Pr), defined in the domain Hr = 0 Br( 0). The conclusion follows from Theorem 

o in the case of (H a ), repeating the proofs of Lemmas 14.614.71 and the arguments of the proof of 
Theorem IPI . 


C. Properties of the global branches In the remaining part of this section, we state some 
further properties of the global branch C\ x , both in the bounded and the unbounded domain 
case. For similar properties possessed by solutions of nondegenerate elliptic equations, we refer to 


Lemma 4.9. Assume that fl is a bounded domain and the condition (Tt a ) is fulfilled, A > Ai is 
fixed and u\, r £ C\ l r and u\ £ C\ 1 . Then, we have that 

u\, r (x) < u\(x), for any x £ fl r , and any r —► 0, 


and 

u\, r * n\ in L^ c (fl \ {0}), osr-»0. 


Proof: The solution u\, r satisfies (P) r , while u\ satisfies 

r —div (cr(x)Vu) = \u — \u \ 21 u, 

l w|an = 0, u\dB r > 0. 

Having in mind, that both u\^ r and u\ are sufficiently smooth and positive functions on f ~l r , from 
the Comparison Principle ins, Theorem 10.5], we conclude the first assertion of Lemma. 

Next, we proceed as in Lemma 14.51 Since u\ tT £ -Dq’ 2 (H, <r), it holds that ||'UA, , r|li, 2 £(n) < 
K ||u.\,r||.Di'2(Q cr), where K is given in 12.21 and is independent of r. Hence, from Lemma |4.1I we 
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have that u\ <r is uniformly bounded in L°°(fl). Consider ip = u — u\ <r . Then, from (391 Theorem 
8 .8] we obtain that 

H^llwfoclsAto}) - c Wi’Ww^(n) +°{r), as r -> 0, 

for some positive constant c independent from r. Then, by a standard bootstrap argument, we 
conclude that u\ tT —> u\ in \ {0}) and the proof is completed, o 

Similar results may be obtained for the unbounded domain case. 

Lemma 4.10. Assume that 12 is an unbounded domain and the condition (Tiff) is fulfilled, A > Ai 
is fixed, u\ ,R e C\ 1R and u\ G C\ 1 . Then, we have that 

u\ ,r(x) < u\(x), for anyx G tt R , 


and 

u\,r->u\ inL“ c (fl\{0}), asR->oo. 

For both, the bounded and the unbounded domain case, we have the following 

Proposition 4.11. Assume that conditions (TL a ) or {Ttpf) hold. Then, 

(i) The global branch C \ 1 bends to the right of Ai (supercritical bifurcation) and it is bounded 
for A bounded. 

(ii) Every solution u G C\ IS is the unique nonnegative solution for the problem mi. 

Proof: (i) Assume that C\ x bends to the left of Ai. Then there exists a pair (A ,u) G R. x 
Z?g’ 2 (fl, cr), 0 < A < Ai, such that 

(4.11) f er(x)|Vu| 2 dx = A f \u \ 2 dx — f |w| 27+2 dx, 

O v O v o 

The last equality implies that 

- A H u lli 2 (n)> with A<Ai, 

which contradicts the variational characterization JUJ of Ai. Thus, Caj must bend to the right of 
Ai. To show that C\ x is bounded for A bounded, we proceed exactly as for the derivation of the 
estimate EZB - 

(ii) Let u G C\ ± , and suppose that v is a nonnegative solution of m with u ^ v. We claim that 
u(x) < v{x), for any x G 12 \ {0}. This is a concequence of Lemma f4.91 since by the Comparison 
Principle we have that 

u\, r (x) (or u\, R (x)) < min{u(x), u(a;)}, 

and of the L)^ c -convergence of u\, r (or u\,r(x)) to u. Then, from 14.1111 we must have that 


f (\u\ 2l v — \v\ 2l u) dx = 0, 
Jn 


which is a contradiction, unless u = v. o 

We emphasize that uniqueness results in the case of semilinear elliptic equations, have been 
treated by many authors. We refer to the discussion in ESI. Theorem 2.4]. For an approach using 
variational methods we refer to P‘21 Theorem 4.1]). 
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5. Convergence to the nonnegative equilibrium, in the case of a bounded domain. 

Theorem ll. ll establish.es for any A > Ai, the existence of a unique nonnegative equilibrium point 
for the semiflow S(t). In the light of Theorem 11.21 in order to prove convergence of solutions of 
(11.11) to the nonnegative equilibrium, it remains to verify (a) that solutions of 111 .Ill remain positive 
for all times and (b) the asymptotic stability of the nonnegative equilibrium. 

Proposition 5.1. Assume that condition ( Tt a ) or (' H < jfi) holds. The set 

T>+ := £ D] ) ' 2 (fl, a) : <f>(x) > 0 on fl| , 

is a positively invariant set for the semiflow S(t). 

Proof: From Proposition roi we have that solutions are globally defined in time. It suffices 
to show that a kind of maximum principle holds, that is, solutions of o> corresponding to 
nonnegative initial data, remain positive. We adapt an argument from Proposition 5.3.1]. Let 
4>o £ -Dq’ 2 (S1, tr), 4> 0 > 0 a.e in ft, and (j) £ C([0, +oo); Z?Q ,2 (fl, a))nC 1 ([0, +oo); L 2 (fi)) the global in 
time solution of o, with initial condition </> 0 . We consider (f> + := max{<^>, 0}, <f> := — min{<^>, 0}. 

Both </>+ and (f~ are nonnegative, (f~, (j) + G^G C([0, +oo); Dl’ 2 (fl, a)) fl C 1 ([0, +oo); L 2 (fl)), and 
we set (j) = (j) + — cf>-. We get from o, that </> satisfies the equation 

(5.1) — div(cr(:r) V0“) — A </>“ + |</>| 27 </>~ = 0. 

Multiplying <EH> by </> and integrating over fl we obtain 

(5-2) [ cr(x)\\/</)-\ 2 dx - A||</>'||| 2 ,+ f H 27 |<r| 2 efe = 0, 

which implies that 

^irii| 2 <ciriii. 

Thus, by Gronwall’s Lemma we obtain 

(5.3) ||0 - (i)|l!a < e ct ||^o ||| 2 = 0, for every t G [0,+oo), 
hence <fi > 0 for all t G (0, +oo), a.e. in fl. o 

Lemma 5.2. Let condition (fH a ) be fulfilled. The unique nonnegative equilibrium point which 
exists for A > Ai is uniformly asymptotically stable. 

We discuss first the stability properties of the zero solution. The linearization about the zero 
solution which is an equilibrium point for any A is 

d t ip — div(cr(a:)Vi/>) — Xi/j = 0, x G fl, 

'f’lon = 0. 

It follows from El, that <fi = 0 is asymptotically stable in Dq’ (fi, a) if A < Ai, and unstable in 
Hg’ 2 (n, cr) if A > Ai. 

The linearization around the nonnegative equilibrium point u of o, is given by 

(5.4) — div(u(a:)V^) — Xip + (2y + l)\u\ 2l ip = 0, 

V’lan = 0: 

and we shall see that for the corresponding eigenvalue problem 

(5.5) — div(cr(a;)V'!/') — Xip + (2y + \)\u\ 2l %p = pip, 

ip\an = 0, 


20 


zero is not an eigenvalue. The weak formulation of lETS!) is 


A(if,u) a := / cr(x)V , 0Vu; dx — A / ifudx 

Jn Jn 

(5.6) + (27 + I) f \u\ 2l ifui dx = /i f i/judx, 

Jn Jn 

for every u> £ _Dq’ 2 (S1, a). The symmetric bilinear form A a : ZDq' 2 (ST, a) x Dq' 2 (Q, a) —y R defines a 
Garding form G pg. 366], since 

,1p) > 11^11^1,2(0 ^ - A||t/?||L2 ( o)- 

Hence, Garding’s inequality is satisfied. Then it follows from Lemmas 12.212.31 and G Theorem 
22.G pg. 369-370], that the problem 15.511 has infinitely many eigenvalues of finite multiplicity, and 
if we count the eigenvalues according to their multiplicity, then 

(5.7) — A < Hi < H 2 < • • • , and fij —> 00 as j —> 00 . 

The smallest eigenvalue can be characterized by the minimization problem 

(5.8) m = min A a (if, ip), ip £ D^’ 2 (n, cr), \\ip \\ L 2 = 1 . 

The j -th eigenvalue, can be characterized by the minimum-maximum principle 


(5.9) u, = min ma xA a (ip,ip). 

J MeCj^GM 

where M = {ip £ Dq' 2 (H,ct) : ||'i/’||z ,2 = 1} and Cj denotes the class of all sets M (~l L with L an 
arbitrary j-dimensional linear subspace of Dq' 2 (fl, a). 

By using similar arguments as for the proof of Lemmas 14.214.61 we may see that for 15.51) . the 
(nontrivial) eigenfunction corresponding to the principal eigenvalue p ,1 is nonnegative, i.e ipi > 0 
a.e. on fl. Since pi,ipi satisfy GO) we get by setting u> = u that 

/ a(x)'Vipi Vu dx — A / ip\u dx + (2y + 1) / \u\ 2l ipiu dx = p\ / ip\u. 

J Q, J Q, J J £1 

On the other hand, by setting v = ip\ to the weak formula (I3T91) we get 


/ cr(x)V'0iVw dx — A / ipiudx + / \u\ 2 ^ip\udx = 0. 

J J £1 J 

Subtracting these equations, we obtain that 

(5.10) 2y I \u\ 2ry uip\dx = /ii / uipidx , 

J Q J Q 

which implies that pi > 0 . o 

Proof of Corollary roi The positivity property of Proposition 15.II and Theorem o imply 
that the solution <p(-,t) converges towards the set of nonnegative solutions of 11 . 21 ) as t —> 00, in 
Dl’ 2 (Q, a). In fact, it follows from Lemma f.5.21 that in the case of A > Ai, for any nonnegative 
initial condition <f> 0 , u>((po) = {it}. On the other hand it is not hard to check, by following the 
computations leading to mu rm . that in the case A < Ai, dist(<S(t)B, {0}) —> 0 as t — > 00, for 
every bounded set B C Dl' 2 (Q, a). In this case, the global attractor A is reduced to {0}. o 


Remark 5.3. (Minimization of the Lyapunov Function) We expect naturally, that the nonnegative 
equilibrium points, minimize the Lyapunov function m 0 a ssume that condition (TL a ) holds. 
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It is not hard to check that J is a bounded from below functional on D(' 2 (Q, a). Assume further 
that A < Ai. Then, the variational characterization of Ai 10 ]) , implies that 


J(<t>) > 


1 

27 + 2 


m\ 


27+2 
l 2 ~<+ 2 ( n)> 


for every </> £ Dq' 2 (fl, a). Hence the trivial solution is the global minimizer of the functional J. 
However, for A > Ai the origin is no longer the global minimizer of the functional: consider the 
function tu\, where u\ is the normalized nonnegative eigenfunction associated to Ai and t > 0 is 
small enough. Then from \d. r l\), we have that 


J(tui) = 


ylMI^ 

t 2 1+ 2 

+ 


Ai t 2 




27 + 2 
(A — Ai) t 2 


27+2 

L 2 -v+ 2 (fi) 


ll u i|li 2 (n) 


^27+2 


(A — Ai) t 2 


l u i|li 2 (n) 


ll w i|li 2 (n) + 27 7- 2 H Ul lli^^+ 2 (fJ) 


< 0 . 


The justification of the Palais-Smale condition follows from Lemma (ii). Then, Ekeland's 
variational principle implies the existence of nontrivial minimizers for J. These minimizers are 
the solutions which belong to the branch C\ ± . 

Actually C\^ is a pitchfork bifurcation of supercritical type, where Principle of Exchange of 
Stability holds. 


Remark 5.4. (Degeneracy exponent) Condition lff.il) can be written as a restriction on the “de¬ 
generacy’'' exponent 


(5.11) 


0 < a < 


2 (1 -7(jV_2)) 

27 +I 


:= a 


This is a restriction on the ‘Tate” of decrease of the diffusion coefficient a near every point z £ 
cr 1 {0}. Unique (since the nonlinearity defines a Lipschitz map) and global in time solutions of 
EJ exist, converging towards a global attractor, if cr{x) decreases more slowly than \x — z\ s , 
6 £ (0, a*], near every point z £ cr 1 {0}. As an example we mention the case N = 2 and 7 = 1 
(cubic nonlinearity) where a* = 2/3. However, as it follows from the discussion in JT3], 15.111) 
(or ra; does not possibly define a critical exponent, concerning the existence of global attractor. 
In the case a > a*, the dynamics related to m. could be investigated through the theory of 
generalised semiflows H21 ESI- As in the case of the damped semilinear wave equation examined 
in G3, where uniqueness of solutions is not assumed, one could possibly prove the existence of a 
global attractor in the case a > a*, under the hypothesis that weak solutions of El satisfy the 
corresponding energy equation. 

Possibly, an interesting issue could be, the extension of the bifuraction result and convergence 
to equilibrium, to the complex evolution equation discussed in ESI Writing the stationary problem 
as a real system for the real part <j>\ and imaginary part <f> 2 , one could observe, that at least in 
the case of the Ginzburg-Landau equation with real coefficients, the stationary problem defines a 
potential system S3, which admits a positive principal eigenvalue Ai. The eigenvalue Ai could 
be a bifurcation point, from which two global branches bifurcate. These branches could consist of 
semitrivial solutions (i.e. solutions of the form ( 4 >i, 0 ) or (O,^))- 


Remark 5.5. (Lack of compactess) Our approach concerning convergence to the equilibrium, which 
combines the charachterization of the global attractor of and global bifurcation theory depends heav¬ 
ily on (H a ) and (Tiff), ensuring compactness of the linear and nonlinear operators involved to our 
study, either in the bounded or unbounded domain case. Thus it is natural to ask if a relaxation 
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of the aforementioned conditions which may give rise to noncompatness, could allow for a gener¬ 
alization of the results of Sections 3-5. 

A starting point, could be the generalization of the result concerning the existence of the global 
attractor. One could assume conditions ( H a ) and {Tiff) for some a £ [0, +oo). As it is noted in 
p~81 . Remark 2.1], if cr £ L) oc (f}), Q C R w , N > 2, satisfies (' H a ) then it satisfies (Hp) for any 
(3 > a and if (7i“) holds, then {Tif 1 ) is valid for any 5 £ [0, a]. 

For example, in the case where S2 = one can also consider as an energy space, the space 
Hq(WL n ,a), defined as the closure of C(f(M. N ) with respect to the norm 

1101 Iff ,<7 = [ cr(x)\V(/)\ 2 dx+ [ |0| 2 dx. 

Jr n Jr n 

The embedding Hq(K n ,o) C L 2 (WL n ) although obviously continuous, is not compact, in the case 
where ( Ti a ) and {Tiff), hold for some a £ (0, 2]. Recall also from Remark \2.4\ that Dq' 2 (M. n , a) is 
not compactly embedded in L 2 (R n ), if a grows less or equal than quadratically at infinity, even in 
the case where ( TL a ) is satisfied for some a £ (0, 2). 

For the existence of global attractors for reaction diffusion equations in unbounded domains, 
representative references include [E : , for semilinear and degenerate (porous medium) parabolic 
equations considered on weighted Sobolev spaces) and m- The latter provides with an effective 
remedy for the lack of compactness of the Sobolev imbeddings, with respect to the existence of 
the global attractor for partial differential equations considered in unbounded domains and in the 
natural phase space. The idea of m is based on the approximation of M. N by a bounded domain 
and on the derivation of suitable estimates for the approximation error of the norm of solutions, 
showing that this approximation error is arbitrary small. These estimates allow for the application 
of the method developed in El which makes use of the energy functional associated to the evolution 
equation (in m a reaction diffusion equation): Asymptotic compactness is shown by passing to 
the limit of the nonlinear term of the energy functional as the error tends to zero, establishing the 
existence of a global atractor in L 2 (WL n ). 

It would be possibly interesting to attempt to apply this method, to the degenerate equation of 
the form ED and investigate if new restrictions could arise between degeneracy, nonlinearity and 
the parameters involved, through the process of the derivation of an appropriate energy functional, 
and the estimation of the relevant estimation errors of the generalised Sobolev norms. 

On the other hand, as it is already mentioned in the introduction with respect to the convergence 
to equilibrium, in the unbounded domain case, one has to deal in general, not only with the lack 
of compactness but also with the possible appearance of infinite distinct translates of a unique 
rest point. Thus, it could be also interesting to investigate if the analysis of E2, could provide a 
framework for the generalization of the convergence result for ED, in the unbounded domain case. 

6. Applications of the Global Bifurcation Result to general elliptic equations 

We conclude, by mentioning some other examples of degenerate elliptic equations for which, 
extensions of the results of Section 4, could be investigated. 

A. Semilinear Equations We consider the semilinear problem 

-V(a{x)Vu) = \f(x)u-g(\,x,u), 

(6.1) 

w|on = 0. 

and the corresponding linear eigenvalue problem 

—V(ct(x)Vm) = A f(x)u, 
u\on = 0, 
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( 6 . 2 ) 


where fl C M. N , N > 2, is an arbitrary domain (bounded or unbounded). In this case, the coefficient 
functions satisfy: 

2 * 

(T) a f is a smooth function, at least (7^(11), for some £ G (0,1), such that / € L 2 ^~ 2 (f2), and 
there exists flf C f2, |ffj£ > 0, such that /(a;) > 0 for all x G 

(£/)“ g is a Caratheodory function, i.e., g(-,x,-) is a continuous function for a.e. x G f2 and 

g(A, -, it) is measurable for all (A, it) G R 2 . Moreover, there exist nonnegative functions c(A) G C(R) 
2 * 

and p(x) G L°°(tt) (f2), such that |<?(A, x, it)| < c(A) p{x) |it| 7 , for all (A, it) G R 2 and almost 

every x G R w . 

Depending on the particular properties of the coefficient functions, the properties of the global 
branch could be represented by those of the corresponding approximating problems. For some 
applications, we refer to El El SOI ESI 

B. Quasilinear Equations We consider quasilinear degenerate elliptic equations of the form 

(6.3) —V(a(x) |Vit| p_2 Vu) = A b(x) |u| p ~ 2 it + f(x) \u\ 1 ~ 1 u 
and 

(6.4) — V(a(x, it) |Vu| p_2 Vit) = A b(x) |u| p_2 u + f(x) |u| 7_1 it, 

where flia a bounded domain of R w , N > 2 and 1 < p < N. A possible treatment could be based 
on the results of 129 ■ Assume that there exists a function v(x) > 0 in fl satisfying 

(A/i) v{x) = 0 or v(x) = oo in a finite subset ZcO, 

(A4) v G Lj oc (Cl), G Lj oc (n) and v~ s G I^Q), for some s G |max{ A, ^-j, oo ^ . 


The coefficient functions a, b , / satisfy: 

{A) a is a smooth function at least C^, for some £ G (0,1) a.e. in fi, such that 


v(x) 


< a(x) < c p(x), for some c > 0. 


( B) b is a nonnegative and smooth function, at least for some £ G (0,1), such that 

beL°°(ty, 

p* 

(J 7 ) f is a smooth function, at least C^(fi), for some £ G (0,1), such that / G i p » _h+1) , where 
Pi ■= N (s+i)-ps and p < 7 + 1 < p* s . 

We consider the weighted Sobolev space Z?Q’ p (fl, u) endowed with the norm 


\ 1 /P 

v(x) |Vu| P j < 00. 

The space Z?Q ,p (fi, u) is a reflexive Banach space, enjoying the following embeddings: 

i) -Dg’ p (fi, v) ^ L p ‘ (fl) continuously for 1 < p* < N, 

ii) Dq’ p ( Cl, v) c —>• L r (Vl) compactly for any r G [l,p*). 

For further properties of these spaces, we refer to ,'29, as well as for the proof of the following 
results. 


11 ^( 0 ^) 
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Lemma 6.1. Assume that conditions (A/2), (A), ( B ) and ( T) hold. Then, the corresponding to 
\ 6 .A) eigenvalue problem 

(6.5) — V(a(a;) |Vit| p_2 Vzt) = A b{x) |u| p_2 u, 

admits a positive principal eigenvalue X\, given by 

Ai = inf [ a{x ) |V</| P dx. 

f n b (x)\ 4 >\ p dx=iJn 

Moreover, Ai is simple with a nonnegative associated eigenfunction ui. In addition, Ai is the only 
eigenvalue with nonnegative associated eigenfunction. 

Lemma 6.2. Assume that conditions (A/2), (.A), ( B) and {IF) hold. Then, any weak solution 
u e v) of ifOD belongs to L°°(fl). Moreover, u £ for some ( £ (0,1), a.e. in ft. 

Proposition 6.3. Assume that conditions (A/2), (A), ( B) and {IF) hold. Then the principal eigen¬ 
value Ai of 16'. 51) . is a bifurcation point of the problem 16'. 51) . 

Pronosition 16.31 could be extended to a global bifurcation result as follows: Assuming that v sat¬ 
isfies in adition, condition (A/i), the principal eigenvalue Ai is a bifurcating point of a global branch. 
We may adapt the same procedure described in Sections 2-4, by considering similar approximating 
problems. It is interesting to note that in this case, Picone’s identity is still applicable. 

Theorem 6.4. Assume that the conditions (A/i), (A/2), (A), ( B ) and {IF) hold. Then the branch 
bifurcating from the principal eigenvalue Ai of 16'. 51) . is a global branch of solutions of the problem 
16'. 51) . Moreover, any solution which belongs in this branch, is nonnegative. 

Concerning problem m, we assume that a(x,u), b{x,u ), are sufficiently smooth functions 
satisfying the following conditions: 

(AtS) < a(x, s) < c g(|s|) v{x), 0 < b{x, s ) < b{x) and lim s _>o o,{x, s) = a(x), lim s ^o b(x, s) = 
b(x) uniformly for a.e. x £ 0. 

Here c > 0, g is a nondecreasing bounded function and a, b satisfy conditions (A) and (£>), 
respectively. 

Based again on J2|j, and the analysis of Sections 2-4, we may prove 

Theorem 6.5. Assume that condition (A<S) holds. Then the principal eigenvalue Ai of 1 6'. 51) . is 
a bifurcation point of the problem UH\ )- Moreover, the corresponding branch is global, and any 
solution which belongs to this branch, is nonnegative. 

The quasilinear problems, 

—V((fc \x\ a + m la;^) |VM| p ~ 2 Vri) = A \u \ p ~ 2 u + f{x) |u| 7-1 m 

and 

—V((k \x\ a + to \x\P) (1 + e _1 /“ )|Vm| p_ 2 Vm) = Ae _ “ \u\ p ~ 2 u + /( x) ImI 7-1 ^, 
for some 0 < a < rninjp, N{p — 1)}, — N < /3 < 0 and k, to. nonnegative constants, could serve as 
examples. 
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